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ABSTRACT 

We present an analysis of the Minkowski Functionals (MFs) describing the WMAP 
three-year temperature maps to place limits on possible levels of primordial non- 
Gaussianity. In particular, we apply perturbative formulae for the MFs to give con- 
straints on the usual non-linear coupling constant /nl- The theoretical predictions are 
found to agree with the MFs of simulated CMB maps including the full effects of radia- 
tive transfer. The agreement is also very good even when the simulation maps include 
various observational artifacts, including the pixel window function, beam smearing, 
inhomogeneous noise and the survey mask. We find accordingly that these analytical 
formulae can be applied directly to observational measurements of /nl without rely- 
ing on non-Gaussian simulations. Considering the bin-to-bin covariance of the MFs in 
WMAP in a chi-square analysis, we find that the primordial non-Gaussianity param- 
eter is constrained to lie in the range —70 < /nl < 91 (95% C.L.) using the Q+V+W 
co- added maps. 

Key words: Cosmology: early Universe - cosmic microwave background - methods: 
statistical - analytical 



1 INTRODUCTION 

The existence of non-Gaussianity in primordial density fields 
has the potential to provide a unique observational probe 
that will enable discrimination among wide variety of in- 
flationary models of the early Universe. Versions of the in- 
flation scenario based on the idea of a single slow-rolling 
scalar field predict levels of non-Gaussianity too small to 
be observed. On the other hand, multi- field inflation models 
and models with a non-standard kinetic term for the infla- 
ton may yield larger non-Gaussian effects which could in 



tions to the curvature perturb ation during the matter era 
(e.g. iKomatsu fc Spergel|[200ll ): 



<£> = ^+/nl(0 2 -(0 2 }), 



(1) 



vations (e. 


e. iBartolo et all |2002|; iBernardeau & Uzan 2002: 


Lvth et al. 


20031: iDvali et all 12004 lArkami-Hamed et al. 


20041: [Alishahiha et al.ll2004l:lBartolo et al.ll2004l:IChen et al. 


200(1: iBattefeld & Battefeldll2007l; Kovama et alj|2007h 



In this paper we focus on the local parametrisation of 
primordial non-Gaussianity by including quadratic correc- 
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where <j> represents an auxiliary random-Gaussian field and 
/nl characterizes the amplitude of the non-linear contribu- 
tion to the overall perturbation. This local form is motivated 
by the simple slow-rolling single scalar inflation scenario and 
other models, including curv aton models; for an alter native 
parameterization of /nl, see (|Creminelli et a l. 2007a). Cur- 
rent observations are not sufficiently sensitive to detect the 
wavelength dependence of /nl so a constant /nl provides a 
reazonable parameterization of the level of non-Gaussianity. 

Analysis of the angular bispectrum for the WMAP 
3-year data provides a constraint o n /nl to lie b e tween 
-54 and 114 at the 95% C.L. (IKomatsu et al l 120031 ; 
ISpergel et al. l l2007l h lCreminelli et al.l (|2007ah obtains more 
string ent constraint —3 6 < / nl < 100. On the other 
hand, lYadav fc Wan dclt ( 20Q3) recently reported a detec- 
tion of primordial non-Gaussianity at greater than 99.5% 
significance. Further detailed analyses of non-Gaussianity is 
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clearly necessary in order to reconcile and understand the 
various constraints and claimed detections. 

Different approaches to the study of non-Gaussianity 
exploit different statistical properties and will be sensitive 
to different aspects of the behaviour of the pattern being 
tested. In general, there is no unique statistic to describe the 
non-Gaussian nature of a sample in a complete manner. A 
given method may have strong discriminatory power for one 
particular form of non-Gaussianity, but this is not necessar- 
ily the case for all possible alternative distributions. Testing 
non-Gaussianity therefore requires a battery of complemen- 
tary techniques rather than a single approach. It is particu- 
larly important to use different statistical approaches in the 
context of primordial non-Gaussianity, because the physical 
mechanism responsible remains unknown. Furthermore, in 
the real world, issues including survey masks and inhomo- 
geneous noise have to be taken into consideration. Different 
statistics may be sensitive to different systematics and fore- 
ground artefacts, so that complementary analysis using dif- 
ferent statistics are essential for a robust detection. Analyses 
using different statistical methods are useful to validate or 
refute basic theoretical models and constrain model parame- 
ters more accurately. The most commonly used statistic for 
non-Gaussian analysis is the bispectrum (or even trispec- 
trum) which focuses on information contained within three- 
point (or four-point) correlations. Other approaches repre- 
sented by Minkowski Functionals (MFs) and genus statistics 
(one of MFs) utilize information concerning the integrated 
morphology and topology of the density structure, and are 
dependent on all order of correlation functions. Their robust- 
ness and generality therefore makes them ideal complements 
to standard correlation analyses. 

In this analysis, we focus on the local model of primor- 
dial non-Gaussianity cha racterized by /nl in Equation (flj. 
ICreminelli et al.l (|2007rJ ) show that the bispectrum is the 
optimal statistics in the estimation of /nl and then other 
statistics (e.g., trispectrum) are useless even if there are dif- 
ferent foreground contaminations. This is, however, only the 
case when the local model exactly describes the real uni- 
verse. Other forms of non-Gaussianity, different from the 
one purely characterized by /nl, which may exist in a real 
observation, could make the fit of the theoretical estima- 
tion as a function of /nl to the observation worse and also 
influence the estimation of /nl among different statistics. 
Different statistical approaches are, therefore, still useful to 
test the assumed model of primordial non-Gaussianity and 
also possible observational systematics by checking if they 
have a reasonable goodness of the fit to observations and 
thus give a consistent limit on /nl compared to that that 
from the bispectrum. 

In this paper, we present a measurement of primor- 
dial non-Gaussianity from the MFs of the WMAP three- 
year temperature m aps. We appl y pert urbative formulae 
recently derived by iHikage et al.l (|2006h to do the com- 
parison with observati ons; previous analyse s rely on non- 
Gaus sian simulations l|Komatsu et al.l 120031 ; ISpergel et al.l 
2007). The agreement of the theoretical predictions with 
non- Gaussian sim ulations has alread y been established in 
Sachs- Wolfe limit (|Hikage et al.|[200r3) . In this paper, we ap- 
ply non-Gaussian simulations based on full radiative trans- 
fer computations and then demonstrate that the analyti- 
cal predictions accurately reproduce the simulation results. 



iGott et all i|2007h already derived an analytical formula for 
the genus statistic in the Sachs- Wolfe approximation to com- 
pare with WMAP data. Our analysis takes more detailed 
physics into account and is consequently more accurately 
applicable to a wider range of scales. 

Observational effects (including antenna beam pattern, 
inhomogeneous noise and the survey mask) could be other 
sources of confusion. From a comparison with simulations 
including these observational issues, we find that the obser- 
vational systematics are negligible to estimate the primor- 
dial non-Gaussianity from WMAP data directly using our 
method. 

The organization of this paper is as follows. In Sj2] the 
WMAP three-year data studied here are briefly introduced. 
In |3]we test whether the perturbative formulae well describe 
the MFs for the non-Gaussian simulation maps even includ- 
ing the various observational effects mentioned above. In Sj4l 
we show the MFs for WMAP three-year temperature maps 
compared with theoretical formulae and give constraints on 
/nl- § [5] is devoted to a summary and the conclusions. 



2 WMAP THREE- YEAR DATA 

The CMB temperature maps derived from the WMAP ob- 
servation are pixelized in HEALPi x format with the total 
number of pixels n p i x = 12iV side dGorski et alj 12005^ . In 
our analysis, we use the maps for Q, V and W frequency 
bands with iV S id e = 512. The linearly co-added maps are 
constructed using an inverse weight of the pixel-noise vari- 
ance o"o /N bs , where cto denotes the pixel noise for each dif- 
ferential assembly (DA) given in lBennett et alj (|2003bft and 
Aobs represents the full-sky average of the effective number 
of observations per each pixel. We adopt two maps with dif- 
ferent combinations of frequency bands: V and W (written as 
"V+W") and Q, V and W (written as "Q+V+W"). The co- 
added maps are masked with the KpO galaxy mas k includ - 
ing point-source mask provided by iBennett et af] (l2003oh . 
which leaves 76.8% of the sky available for the data analy- 
sis. 

In comparison with WMAP observations to give con- 
straint on /nl in |4] a ACDM cosmology is assumed with the 
cosmological parameters at the maximum li kelihood peak 
from the WMAP three-year data only fit (|Spergel et al.l 
12007ft : Q b = 0.04309, fi c dm = 0.211, fl A = 0.74591, H = 
71.227 km s" 1 Mpc" 1 , r = 0.08982, and n s = 0.95537. The 
amplitude of the primordial fluctuations has been normal- 
ized by the first acoustic peak of the temperature power 
spectrum, 1(1 + 1)C i/(2-k) = 5617.05(/xK) 2 at I = 220 
l|Hinshaw et al.ll2007ft . 



3 PERTURBATIVE FORMULAE VERSUS 
NON-GAUSSIAN SIMULATIONS 

3.1 Perturbative Formulae of MFs for CMB with 
Primordial Non-Gaussianity 

The topology of random fluctuation fields is generally stud- 
ied using their excursion sets, i.e. regions where the field 
exceeds some threshold level. In a two-dimensional ran- 
dom field such as a CMB temperature map, three MFs 
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are defined: the fraction of area Vb exceeding the thresh- 
old, the total circumference Vi of all the entire exc ursion 
set, a nd the corresponding Euler Characteristic V2 (|Colesl 
1988). We measure MFs for CMB temperature maps as a 
function of the threshold density u, defined as the tempera- 
ture fluctuation AT /T normalized by its standard deviation 
(To = {(AT/T) 2 ) 1 / 2 . Based on the general formalism of per- 
turbat ion theory for MFs |Matsubarall2003h . Ifflkage et ail 
(2006) derived perturbative formulae of the MFs as a func- 
tion of the non- linear coupling parameter /nl (eq. [T]). 

The MFs are separately written with the amplitude and 
the function of v as follows. 



V k (y) = A k v k (u). 



The amplitude Ak, which is determined only by the angular 
power spectrum Ci, is given by 



( 2w )(fc+D/2 u 2 - k u k \y/2a 
1 x - 



4tt L^i 



(21 + 1) [l{l + i)Y awf 



(3) 
(4) 



where uo k = % k ^ 2 /Y{k/2 + 1) gives wo = 1, wi = 2, L02 = 
7r and Wi represents the smoothing kernel determined by 
the pixel and beam window functions and any additional 
smoothing (e.g. a Gaussian kernel). In weakly non-Gaussian 
fields, the function v k (v) can be divided into the Gaussian 
term v k and the non-Gaussian term at lowest order Av k : 



v k (v)=v k >(v) + Av k (v, /nl). 
Each term has the following form 

(G) -v 2 /2tt 1 \ 



(5) 



(6) 



Av k (u, /nl) 



+ 



e- 2 / 2 { 

fc(fc - 1) g (2) 



(i 



(to + OK 2 ) I , (7) 



where H n (v) represent the n-th Hermite polynomials and 
the sk ewness parameters are given in Equations [27- 
29] of iHikage etafl l|2006l) . The amplitude A fc (eq. [3]) is 
not directly relevant to non-Gaussianity but is dependent 
on the shape of Ci. We therefore concentrate on the non- 
Gaussian term Av k hereafter. The quantity Av k is the same 
as the relative d i fferen ce of MFs, which are plotted in Fig. 2 
in IHikage et al.l (|2006l ). except for its normalization factor; 
in this paper the difference of MFs is normalized by A k ([3]) , 
w hile the maximum va lue of MFs for Gaussian fields is used 
in IHikage et~aH i|2006h . 

3.2 Comparison with Non-Gaussian Simulations 

The above analytical formulae have already been found to 
match accurately the MFs for non-Gaussian m aps in Sachs- 
Wolfe limit (Appendix C in lHikage et al.ll2006l ). Here we test 
them against non- Gaussia n simulations including th e full ra- 
diative transfer function (ILiguori et ali i 2003, 2007J). As we 
mentioned in the introduction, actual observations of CMB 
also involve different effects which may produce other confu- 
sions: the pixel window function, beam smearing, the inho- 
mogeneous noise, survey mask and so on. We include these 



observational effects into the simulations to check whether 
they could have a systematic effect on our topological mea- 
sures. 

The cosmology in the non-Gaussian simulations is 
based on Lambda CDM, but the cosmological param- 
eters have slightly different values from WMAP three- 
year best-fit; Q b — 0.05, f2 dm = 0.25, Qa = 0.7, Ho = 
65 km s _1 Mpc -1 , r = 0, and n s — 1. The amplitude of pri- 
mordial fluctuations is, however, set to be same as WMAP 
three-year best-fit value 1(1 + l)C l /(2-n) = 5617.05(^K) 2 . 

Observational effects related to WMAP data are in- 
cluded as follows. First we convolve the original simulation 
maps with the Q+V+W co-added beam transfer function 
with inverse weight of the full-sky averaged pixel-noise vari- 
ance in each DA. Next we add independent Gaussian noise 
realizations following the noise pattern co-added with the 
same weight. The simulation map is then masked with the 
KpO Galaxy mask. Finally we smooth the simulation maps 
using a Gaussian filter with a smoothing scale of 8 S , 



Wi = exp 



-\i(i + i)el 



(8) 



The MFs are sensitive to the resolution (or smoothing) 
scale of a density field and thereby we can obtain a vari- 
ety of information from density fields using different lev- 
els of smoothing. The information extracted from varying 
smoothing scales is nevertheless limited because they are all 
derived from the same original field; the smoothed fields are 
not independent. Here we focus on the field smoothed by 
three different smoothing scales 10', 20' and 40', where the 
limit on /nl is sufficiently converged. To remove the effect 
of the survey mask near the boundary of the mask, we only 
use the pixels more than 29 s away from the boundary. The 
sky fraction used in the analysis for each smoothing scale is 
41% for 9s = 40', 62% for 8 3 = 20' and 73% for 6 S = 10'. 

The MFs for the measured CMB temperature 
anisotropy are computed from the integral of the curvature 
of iso-temperature c ontour lengths (the details are described 
in Appendix A.l. of Hik age et al.l l2006). The binning range 
of v is set to be —3.6 to 3.6 with 18 equally spaced bins of 
v per each MF. This binning way produces well converged 
results irrespective of other choices of the range of v and the 
number of bins. 

We obtain the normalized MFs (eq. [3]) with the am- 
plitude A k computed from Ci of each realization. Then the 
residuals of the normalized MFs from Gaussian predictions 
Av a are calculated at each bin of a, which denote a thresh- 
old value v, a kind of MF k, and a smoothing scale parame- 
terized with 6 S or iVside. Even when the MFs of Gaussian re- 
alizations are computed, however, the function Av a are not 
exactly equal to due to the effect of pixelization, survey 
mask and other numerical artifacts. We therefore measure 
the deviations from the average of the measurements over 
Gaussian realizations and subtract them as 



Av a — Av a — [Av a ] Gaussian, mean ■ 



(9) 



In Fig. [T] we compare the analytical predictions of vari- 
ance, skewness, and MFs with the measurements from the 
simulations for /nl = 100. The simulated results are the av- 
erage over 200 realizations and the error-bars represent the 
error for the average (the sample variance divided by the 
square-root-of 200). The averaged measurements for Gaus- 
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sian CMB maps are subtracted from those for non-Gaussian 
maps in the simulated plots including variance and skewness 
as well as MFs (see eq. [9]). In the plots, we adopt the Gaus- 
sian maps which are generated from the same realizations 
of linear potential fields (0 in eq. [1]) as the non-Gaussian 
CMB maps. The sample variances, represented by the error- 
bars, are cancelled very well in such plots so one can focus 
on the systematic effect of primordial enon-Gaussianity. The 
analytical formulae are found to agree with the simulations 
extremely well even including all observational effects. This 
indicates that we can measure /nl from direct comparison 
of the analytical formulae with observations without hav- 
ing to worry excessively about the presence of such sys- 
tematics. We also check that both the artificial systematics 
[Au Q ] Gaussian, mean and covariance matrix are not strongly 
dependent on the details of cosmology. These results are en- 
couraging, but not unexpected: being based on integrated 
properties, the Minkowski Functionals are expected to be 
robust to such effects. 



4 CONSTRAINTS ON PRIMORDIAL 
NON-GAUSSIANITY FROM WMAP 
THREE- YEAR DATA 

4.1 Covariance Matrix for MFs 

We have adopted a maximum likelihood method to esti- 
mate the best-fit value of /nl and its associated uncer- 
tainty. In 'nearly Gaussian' fields, the distribution functions 
of Av a are well described as multivariate Gaussians. The 
likelihood function of /nl is, therefore, simply proportional 
to exp(— x 2 (/nl)/2) where X 2 (/nl) is computed using the 
theoretical formulae (eq. [7]) as 

X 2 (/nl) = ^[A^-A«£ thcory '(/ N L)]E^, 



distribution of the true value /^'™ c ' . 



Here we consider that 



x [A«i° bs) -A^ ncow (/ N L)] ; 



(theory) / 



(10) 



where a and a' denote the binning number of threshold val- 
ues v, different kinds of MF k, and smoothing scale param- 
eterized with 8 3 or iV s id e - The full covariance matrix T, aa i is 
required because MFs are strongly correlated between dif- 
ferent v, different kinds of MF, and also different ./Vside or 6 S - 
We estimate the covariance matrix of MFs from 1000 Gaus- 
sian simulation maps including the pixel and beam window 
function, KpO survey mask, and inhomogeneous noise for 
WMAP three-year maps. 

The MFs contain information about fluctuations at dif- 
ferent scales, so the results depend on the choice of window 
function. Here, two different types of window functions are 
adopted (in addition to the beam window functions). One is 
a Gaussian window function with the scale characterized by 
8 3 (which is chosen to be sufficiently large compared with 
the pixel size). The other is just the pixel window function 
in HEALPix format with a scale characterized by iVside. The 
multipole components with I higher than 27Vside are cut be- 
cause they suffer from serious aliasing effects. 

Before applying these ideas to the observational data, 
we check if our method based on x 2 analysis is valid using 
simulations with primordial non-Gaussianity. We apply the 
non- Gaussian simulations to check that the likelihood func- 
tion using the equation (|10fl reproduces a valid probability 



the likelihood function of / 



(true) 
NL 



in each realization follows 



a Gaussian distribution around the best-fit value / 



(best) 



p, f (true) , ,(bcstK _ 

Wnl Unl ) ~ 



°7n 



27rCr /N 



■ exp 



1 / j-(true) 
1 1 JNL 



/•(best) 
JNL 



dv a 

. a. a' 



-1 -1/2 



dfs 



where the binning number a represents a threshold v for fc-th 
MF at a given scale iVside (or 6 S if the Gaussian smoothing is 
added) and the covariance matrix E aa / is numerically esti- 
mated from Gaussian simulations. The function 8V a /dfNL 
is independent of /nl (see equation [7]) and thus the un- 
certainty <t/ nl is independent of /nl- According to Bayes' 
theorem, /jjl should distribute around /^' L °°' in the same 
way: 



(best) I /-(true) \ 

I/nl ) 



(true) I /"(best) \ 
NL I /NL J 



(13) 



We estimate the distribution function of /|, b L cst ' from 200 



non- Gaussian CMB simulated maps at a given /, 



(true) 
NL 



and 



then compare with the equation (|11[) . The simulated maps 
include observational effects represented by pixel and beam 
window functions, noise, and KpO survey cut for WMAP 
three- year data. Fig. [2] shows the theoretical predictions of 
^(/nT^I/nl^) at flT ] = (solid) and 100 (dotted), 
from the MFs for the combined maps at Gaussian smoothing 
scales S = 20' and 10' where the uncertainty is cr/ NL = 44. 
The histograms show the distribution of /nl from 200 
(non-)Gaussian realizations. The averages of the best fit val- 
ues of /nl from the simulations are respectively ± 44 (for 

/nl UC> = °) and 101 ± 46 ( for /nl UC) = 100 )- The simula- 
tions reproduce the theoretical predictions of the likelihood 
function very well. Our method is thus well established to 
give constraints on /nl from WMAP three-year map. 



4.2 Constraints on /nl from Minkowski 
Functionals for WMAP Three- Year 
Temperature Anisotropy 

The three MFs for the CMB temperature maps from WMAP 
three-year data are respectively plotted with symbols in each 
column of Fig. [3] (left three columns for the "Q+V+W" map 
and right three columns for the "V+W" map). In each col- 
umn, the top panel shows the MF V& at a representative 
scale (6s = 20'). and the lower three panels illustrate Avt 
at different 9 S = 10', 20' and 40'. The perturbative formulae 
with the best-fit value of /nl to each observed MF are plot- 
ted with lines. The best-fit values and la uncertainty are 
written in the left-bottom side of each panel. In top pan- 
els, all of the amplitude of observed MFs are found to be 
smaller than the theoretical estimations. This comes from 
the deficit of the observed power at low I which generates 
the larger amplit ude of MFs determ ined by cri/ero (eq.[3]), 
as pointed out bv lGott et ail (|2007T l. Fig. H shows the same 
plot but for the MFs with the pixel window function only. 
The top panel shows the MF at iV S ide = 128 and the lower 
three panels illustrate Avk for A r s idc=256, 128, and 64. It 
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Figure 1. Comparison between the analytical predictions (lines) and the numerical estimations averaged over 200 realizations of non- 
Gaussian simulation maps (symbols) for /nl = 100; Upper-left: variances ctq and o\ (eq. [4]) Upper-right: skewness parameters S"' a ' 
(a = 0, 1, and 2) in the equation J7J, Middle: MFs for non-Gaussian fields, Vk (eq. [2]), Lower: the difference ratio of MFs Avk 
(eq. [7]). CMB maps are smoothed with a Gaussian kernel Wi = exp[— 1(1 + l)£>^/2] where 9 a denotes the smoothing scale. The fully 
radiative transfer function is considered for both the theoretical predictions and the simulations. The simulations also include the various 
observational effects for WMAP three-year Q+V+W coadded map; pixel window function, beam smearing, inhomogeneous noise pattern, 
and KpO cut. The error-bars represent the errors for the averaged simulation results over 200 realizations (the sample variance divided 
by the square-root-of 200). 
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Figure 2. The distribution function of the best-fit value of /nl using WMAP three-year mock simulation maps (histogram). We use 200 
realizations of Gaussian simulations (solid) and non-Gaussian simulations with /jj™ 6 ' = 100 (dotted) respectively. The best-fit values are 
obtained by fitting the analytical formulae (eq. [7]) to all of the MFs for the simulations at 8 S = 10' and 20' combined. For comparison, 
we plot the likelihood function of /nl at /^'™°' = and 100 with o/ NL = 44 (eq. 1111 ). which is the expected uncertainty of /hl from 
all of the MFs for CMB maps at 6 S = 10' and 20' combined. 



is interesting that all MFs at A r s id e =64 have large positive 
values of /nl, though the significance is less than 2a. 

Table [1] lists the best-fit values and the la uncertainty 
of /nl for each MF and their combined values at differ- 
ent sets of Gaussian smoothing scales 9 S . The la uncer- 
tainty of /nl is estimated from the range of /nl with 
A\ 2 = X 2 ~ Xmin < 1- The minimum of chi-square Xmin 
and the goodness-of-fit P x 2 >x 2 are listed for each fit. The 
results for the pixel window function only are shown in Ta- 
ble [21 The goodness-of-fit values are reasonable for all the 
fits, which means that the simple form of the primordial 
non-Gaussianity (equation [T]) well describes the behaviour 
of the observed MFs. In other words, present observations 
are too uncertain to allow the extraction of any further in- 
formation about primordial non-Gaussianity (e.g. scale de- 
pendence of /nl). The constraint —70 < /nl < 91 at 95% 
C.L. is obtained from all MFs for the Q+V+W co-added 
map at combined different Gaussian smoothing scales of 10, 
20 and 40 arcmin. A similar constraint is obtained from the 
MFs with pixel-window only as —84 < /nl < 105. The re- 
sults from Q+ V+W co-added map are consistent with th e 
previous ones (jSpereel et al ] |2007l : ICreminelli et al-lfcoOTah . 

There is some frict ion (but not disagreemen t) between 
our results and those bv lYadav fc Wandeltl (|2008T l : our V+W 
analysis finds /nl = —22 ± 43 whereas they find /nl = 
87+30. Moreover our averaged /nl decreases from Q+V+W 
to V+W whereas their /nl increases. This is very interest- 
ing because there is the possibility that e.g., foregrounds 
and point sources might be biasing one of the two results. 
lYadav fe Wandeltl l)2008l ) show in their analysis that these 
effects do not seem to contaminate the primordial bispec- 
trum measurement significantly. It will be then important 
to check their effect on the MFs statistics in order to verify 
if this can explain the differencies among the two results. 
However the observed discrepancies show already how ana- 



lyzing non-Gaussianity using different statistics can provide 
additional interesting information. 



5 SUMMARY AND CONCLUSIONS 

We have presented an analysis of MFs for WMAP the 
three-year temperature maps to limit the primordial non- 
Gaussianity characterized by the nonlinear coupling param- 
eter /nl • To do this we compared perturbative formulae for 
MFs of weakly non-Gaussian fields directly with the obser- 
vations. The analytical formulae are found to be in excel- 
lent agreement with results from non-Gaussian simulations 
of CMB maps including full radiative transfer effects. The 
agreement is still very good when including systematic ob- 
servational effects including the KpO survey mask, pixel and 
beam window functions, and inhomogeneous noise distribu- 
tion for WMAP three-year data. 

We have performed a \ 2 analysis to the comparison of 
the analytical formulae with WMAP three-year data. The 
fits of the analytical formulae to the observations are ac- 
ceptable and we thus obtain a robust constraint of —70 < 
/nl < 91 at 95% C.L. from the Q+V+W coadded maps with 
Gaussian filter at different scales 10', 20' and 40' combined. 
The r esult is consistent with previous results (ISpergel et alj 
120071 ; ICreminelli et al.ll2007al ; lYadav fc WandeltboOct ). 

The behaviour of the results for the V+W maps raises 
some interesting issues; o ur constraint is neg a tively shifted 
-108 < /nl < 64 while lYadav fc Wandeltl (|2008m find a 
more positive range 27 < /nl < 147. The difference between 
the two results should be clearer in the near future survey 
represented by Planck. It is worth investigating this result 
in further detail through a careful analysis of foregrounds 
and point source effects. This will be the subject of future 
work. 
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Figure 3. Comparison between MFs for WMAP three-year temperature maps (symbols) and the analytical formulae with the best-fit 
value of /nl f° r each MF (lines). The MFs are calculated from the Q+V+W co-added map (Left) and the V+W map (Right). Top panels 
show the MFs Vk (fe=0,l, and 2) at 8 a =20 arcmin, and other panels illustrate Auj, (eq. [7]) at 6* s =10, 20 and 40 arcmin respectively. 
Error-bars denote the standard deviation of MFs at each bin of v computed from 1000 Gaussian realizations including the WMAP 
three-year noise distribution, KpO mask and pixel and beam window function. The systematics due to the pixelization effect is estimated 
from the Gaussian realizations and is subtracted from the observed MFs (see eq. [9]). 
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Table 1. The constraints on /nl from MFs for WMAP three-year co-added maps, Q+V+W and V+W, at different Gaussian smoothing 
scales 9 S [arcmin] of 40, 20 and 10 and their combination. For the calculation of their constraints, the maximum likelihood method are 
employed with the analytical formulae (eq. [7]) and the covariance matrix of MFs estimated from 1000 Gaussian realizations. The range 
of v is set from —3.6 to 3.6 with the binning number per each MF of 18. The goodness-of-fit of the analytical formulae is represented by 
the minimum values of x 2 value, x^ in , and the probability with x 2 larger than Xmm- 
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Table 2. Same as Table [T] but for different JV s id e = 64, 128 and 256 without Gaussian smoothing 
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